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Abstract 

Standard formulae of classical electromagnetism for the forces between electric 
charges in motion derived from retarded potentials are compared with those ob- 
tained from a recently developed relativistic classical electrodynamic theory with an 
instantaneous inter-charge force. Problems discussed include small angle Ruther- 
ford scattering, Jackson's recent 'torque paradox' and circular Keplerian orbits. 
Results consistent with special relativity are obtained only with an instantaneous 
interaction. The impossiblity of stable circular motion with retarded fields in either 
classical electromagnetism or Newtonian gravitation is demonstrated. 



Keywords; Special Relativity, Classical Electrodynamics. 
PACS 03.50.De 03.30. +p 



1 Introduction 



Classical electromagnetism (CEM), in particular Maxwell's equations, played an im- 
portant role in the development of special relativity (SR). This docs not however mean 
that CEM, even as presented in modern text books, is a theory fully compatible with 
SR and and devoid of interpretational problems. This point was strongly emphasised by 
Feynman [1]. On the other hand, the most successful physical theory of the 20th cen- 
tury, as much for its powers of prediction as for its precise experimental verification, is 
quantum electrodynamics (QED) [2] . Again, as stressed by Feynman [3], the elements 
of this theory are very simple: charged particles and real or virtual photons described 
by appropriate wavcfunctions or propagators and a single coupling constant, the elemen- 
tary electric charge. In QED the mechanical forces between charges are mediated by the 
exchange of space-like virtual photons. As discussed, in detail, in a recent paper by the 
present author [4], the exchange of such photons is instantaneous in the centre-of-mass 
(CM) frame of two interacting charges. This is a consequence of the relativistic formula 
for velocity in terms of kinematical quantities: v = pc^/E. Energy-momentum conserva- 
tion shows that E(7*) vanishes whereas, in general, p(7*) is non-zero for a virtual photon, 
7*, in the CM frame of any scattering process, implying that its velocity is infinite. This 
conclusion is confirmed by the study, in [4], of the invariant amplitude in momentum 
space for M0ller scattering, e'e" — > e~e~, and its Fourier transform. Recently, convincing 
experimental evidence has been obtained [5, 6] 

for the non-retarded nature of 'bound' magnetic fields with r^^ dependence, (associ- 
ated in QED with virtual photon exchange) in a modern version, probing small r values, 
of the Hertz experiment [7] in which the electromagnetic waves associated with the prop- 
agation of real photons (fields with r"^ dependence) were originally discovered. 

There is no reason to suppose that the space-like virtual photons responsible for the 
electromagnetic forces in M0ller scattering should not also be at the origin of the elec- 
tromagnetic forces between the electrons in two current-carrying conductors separated by 
macroscopic distances. Thus there should be a close connection between CEM and QED. 
Just this connection was explored in the second part of [4]. In the first part of this paper 
all of the mechanical equations of CEM, as well as Maxwell's equations, were derived 
from Coulomb's inverse-square force law, SR and Hamilton's Principle. The equations 
describing inter-charge forces, although consistent with standard CEM formulae in the 
limit where 0(/3^) and higher correction terms are neglected, differ significantly from them 
when such terms are included. The main aim of the present paper is a comparison of the 
RCED (standing for Relativistic Classical Electro-Dynamics) formulae from [4] with the 
conventional text book formulae of CEM, to be found, for example, in [8, 9, 10]. In- 
deed it is found that the latter, corresponding, for moving charges, to retarded potentials 
and fields do not respect the constraints of SR. As will be seen, this is one source of 
the previously unresolved problems of CEM which Feynman found to be of such great 
interest. 

The other source of problems and paradoxes in CEM is that electromagnetic forces 
apparently do not respect Newton's third law. An explanation in terms of non-mechanical 
momentum and angular momentum residing in static electric and magnetic fields has been 
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attempted by many authors [11, 12, 13, 14, 15, 16, 17] over the last half-century. The 
present author's opinion is that this is not the correct explanation, and that Newton's 
third law is actually always respected by electromagnetic forces, without the need to 
introduce non-mechanical degrees of freedom. One important point in this connection is 
the instantaneous nature of electromagnetic forces, a second is the correct treatment of 
relativistic effects, and the third is the understanding that the mass of an object is only 
a constant in the complete absence of interactions, which is evidently not the case for 
objects undergoing mutual electromagnetic interactions^ . 

In fact it is clear that when all physical quantities are properly defined Newton's third 
law must be obeyed because it is a necessary consequence of Newton's first law, (the 
law of inertia) and the definition of 'force' provided by the relativistic generalisation of 
Newton's second law. The proof of this is now briefly sketched. Given a dynamical theory 
it is always possible to obtain an equation such as: 

— * 

Po^ fo{X,ai,aj...) (1.1) 

where po is the relativistic momentum (see (1.6) below) of an object O, the dot denotes 
a time derivative and fo is a known function of the space-time position X = {ct;x) of 

— * 

O and some fixed parameters ai,aj.... The vector fo is the relativistic generalisation 
of what is defined by Newton's second law as a 'force'. Actually, however, this naming 
operation is the entire physical and logical content of the second law. The dynamics is 
completely defined by the knowledge of the time derivative of the relativistic momentum 

(1.1) . Nothing is gained, except in an anthropomorphic sense, by baptising the right side 

— * 

of this equation a 'force'. If it happens that /o = in (1.1) then: 

Po — constant (1.2) 
which is Newton's first law. Suppose now O is composed of two parts Oi and O2 so that: 

Po = POi + P02 = constant (1.3) 
Differentiating (1.3) with respect to time and transposing gives: 

Poi = -Po, (1-4) 

or, using (1.1) 

fo, = -fo, (1-5) 

which is Newton's third law. It is clear that it is a necessary consequence of the first law 

(1.2) and the definition (1.1) of 'force' provided by the second law. The third law must 
then be obeyed if the relativistic momentum is correctly defined. 

The general definition of the relativistic momentum of an object is: 

p = m{X)^^ — m{X)^v (1.6) 
or 



^The concept of an effective mass for a particle, depending upon the strength of the interactions it 
undergoes, is a familiar one in the theory of conduction in the solid state. The relation between the 
energy and momentum of a conduction electron can be markedly different from that for a free electron, 
for those occupying energy levels near to the upper or lower boundaries of a conduction band [18]. 
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where v is the object's velocity, r its proper time, 7 = I/a/I — (3 = v/c and X is the 
space-time position of the object. Only if the latter is in free space, and undergoes no 
interaction with the force fields of other objects, is the mass, m, a constant. A simple 
example of the variability of the masses of interacting objects is provided by the case of 
two equally massive, equally charged, objects moving towards each other along a straight 
line in their common CM frame. At a certain distance, xq, of each object from their 
common center of energy, they come to rest, due to their mutual repulsion. The energy 
of the system in the CM system (i.e. its mass times c^) is then: 

E = 2mc^ + ^ = 2m{xo)c^ (1.7) 

where q is the charge of each object and m its mass in the absence of interactions. The 
definition of the spatially dependent mass in the last member of (1.7) follows from the 
symmetry of the problem -a system at rest composed of two equal sub-systems at rest. 
Thus: ^ 

m{xo)-m+^ (1.8) 

For larger values, x > xq, of the distance, x, of each object from the centre of energy, 
■m{x) will vary in such a manner as to respect the conservation of relativistic energy: 

= 4((m + ^) V + pV) = A{m{x)\' + p'^c^) (1.9) 

Since E is constant, due to the conservative nature of the interaction, it is clear that 
any change in the momentum of the objects must be accompanied by corresponding 
changes in their masses. Another example of position-dependent masses is given below in 
Section 6 where circular Keplerian orbits of a system of two electrically charged objects 
are discussed. 

The structure of this paper is as follows: In the following section the RCED formulae 
for the electric and magnetic force fields of a uniformly moving electric charge are derived 
from Coulomb's Law and the 4-vector character of the electromagnetic potential. In 
Section 3 small angle Rutherford scattering is calculated in different inertial frames, using 
either RCED force fields or the conventional retarded fields. It is shown that the latter, 
when used in the initial rest frame of an electron scattered by a heavy charged particle, 
give a result incompatible with both the known Rutherford scattering formula and the 
Bethe-Bloch energy-loss relation. In this section it is also shown that the forces predicted 
by the retarded fields do not, in general, respect the law of conservation of relativistic 
transverse momentum. In Section 4 the breakdown of Gauss' law for the total flux of the 
RCED electric force field of a charge in uniform motion is demonstrated. In Section 5 
the recent 'torque paradox' of Jackson is resolved by employing the RCED force fields, 
that, unlike the conventional retarded ones, arc consistent with SR. Section 6 presents an 
analysis of circular Keplerian orbits in RCED. The relativistic generalisation of Kepler's 
third law for such orbits is derived. In Section 7 the question raised by 19th century 
astronomers, and recalled by Eddington [19], concerning the possibility of stable orbital 
motion under the infiuence of retarded fields, is addressed, for the case of circular orbits 
in electrically bound systems, and a negative answer is found. The final section is a 
summary of the results obtained. 
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2 Force fields of a moving charge in relativist ic clas- 
sical electrodynamics 



The expressions for the electric and magnetic fields in RCED found in [4] are here 
rederived for the simpler case of a uniformly moving charge by exploiting the 4-vector 
character of the electromagnetic potential, A, and assuming Coulomb's law of force be- 
tween static charges with an instantaneous interaction. The usual covariant definitions^ 
of the electric and magnetic fields ^[20] are used: 

— * 

E = -VAo--^ (2.1) 
c at 

B ^ V xA (2.2) 

Denoting by S' (coordinates t', x') the rest frame of an electric charge q and S (coordinates 
x) a frame, moving with velocity -v along the x'^ axis, in which the fields are to be 

evaluated, then A' = (g/r; 0,0,0) in S'. If X, X' are space-time points in S and S' 
respectively at which the fields are to be evaluated while Xq and X'^ are the corresponding 
space-time positions of the charge, the quantity r is the invariant space-like interval given 
in S' by^: 

r2 = -{X' - X'^f (2.3) 
In view of the instantaneous nature of the Coulomb interaction: t' — t'^, so that: 

= (x' - x\f (2.4) 
The Lorentz transformation equations between the frames S' and S: 

A = 7(^'i + M) (2.5) 
^0 = 7(>lo + M) (2.6) 
give, remembering the Lorentz-scalar character of r: 

^ = (^;^.o,o) (2.7) 

r r 

Substituting (2.7) into (2.1) and (2.2) , and performing the partial differentiations, noting 
that, for the partial time derivative in (2.1), x is held constant while is allowed to vary, 
the following expressions for the RCED electric and magnetic force fields are obtained: 

E(RCED) = (2.8) 

B(RCED) = ^^i^' (2.9) 

where j = qu = qljc; jv, 0, 0) and u is the 4-vector velocity of the charge q. Choosing r 
in the 12 plane and with unit vectors i, j parallel to Oxi and 0^2 respectively, (2.8) can 
be written in the convenient form: 

E(RCED,Se) = ^(^^ + 7isin^) (2.10) 



^Gaussian electromagnetic units are employed. 
^ These are derived from first principles in [4] 
time-like metric is used for 4-vector products 
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where cos^ = {f ■ v)/\r\\v\. This is the RCED formula used in the following section. 

3 Small angle Rutherford scattering with instanta- 
neous or retarded fields 



Rutherford scattering of an electron of charge, — e, and mass, m, by a nucleus of 
charge, Ze, and mass, M, is considered in two frames of reference. The first is the rest 
frame, S^r (coordinates (t';x',y',z')) of the nucleus (essentially the overall centre-of-mass 
frame of the two body scattering problem) the second is the rest frame, Sg (coordinates 
{t;x,y,z)), of the initial electron. As in a similar discussion in [21] a large impact pa- 
rameter, 6, is considered so that the geometrical deflection of the electron trajectory in 
Sjv and the change in position of the electron in Sg, during the scattering interaction, are 
neglected. Also neglected are changes in the velocities of the electron or nucleus during 
the scattering interaction. Geometrical parameters in the frames S ^ and S,, are defined in 
Figs, la and lb respectively, while Fig.lc shows the definition of the scattering angle, 9, 
of the electron momentum vector in the centre-of-mass frame Scm — Sn- The scattering 
angle is given, in the small angle approximation that is made here, as : 

e ^ (3.1) 

p 

p is the electron momentum in Scm and Ap is calculated by integrating the transverse 
impulse produced by the electric field of the nucleus during the passage of the electron in 
Sat or of the nucleus in Sg. Small effects due to magnetic fields in Se (which, since they are 
transverse to the electron velocity, do not contribute to |Ap|) are neglected. The electric 
field is calculated either from the RCED formula (2.10) or by use of the standard formula 
for the electric field of a moving charge given in textbooks on classical electrodynamics. 
This is the same as the 'present time' version of the retarded Lienard-Wiechert (LW) [22] 
field, which was first derived by Heaviside [23] . 

In the frame Sat, the electric field at the position of the electron is given, in both cases, 
by a static Coulomb field, since the source charge is at rest: 

E(Sn, RCED) = E(Sn, LW) = i'E^iSN) + j'Et(Sn) = ii' cos i) + / sin i\)) (3.2) 

where is the angle between the vectors r and v (see Fig. la) and j' are unit vectors 
in the x'- and y'-directions. In the frame Sg, where the nucleus is in motion, different 
formulae are found for the electric field in the two cases: 

EiS,,RCED) = ^(i^ + ^jsinV') (3.3) 
wo ^T^^ _ (zcos^ + jsin^) 

^^^-"^"^^ - ^(l-/3^sin^^)* ^'-'^ 
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where ? and j are unit vectors in the x- and y-directions. Although the transverse fields 
given by (3.3) and (3.4) agree when ip = 7t/2, the angular dependence for fixed r is very 
different, as may be seen in Figs. 2 and 3. which show respectively E^r'^/Ze and Exr'^/Ze, 
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where El and Et are the longitudinal and transverse components of E, repectively, as a 
function of ip, for various values oi (3. It can be seen in Fig. 3 that the denominator of the 
right side of (3.4) strongly damps the strength of the transverse electric field, except for 
values of close to 7r/2, for values of j3 close to unity. On performing the integral over ijj 
to calculate Ap a much smaller value of 9 is then to be expected for a given value of b for 

— * — # 

E{Se, LW) as compared to E{Se, RCED). This is indeed found to be the case. 




Figure 1: Geometrical parameters for the calculation of Rutherford scattering in different 
frames. In a) the nucleus, N, is at rest in Sn, and the scattered electron, e, is in motion. 
In b) the electron is initially at rest in Se while the nucleus is in motion, c) shows the 
dehnition of the electron scattering angle, 6, in the overall center-of-mass frame, Scm- To 
a good approximation, Sn and Sqm sure the same. 
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Figure 2: The longitudinal electric field of a uniformly moving electric charge, scaled by 
the factor r'^/Ze, as a function of angular position ijj, for different relativistic velocities 
(3 = v/c. The slowly varying cosine curves correspond to the RCED formula (3.3), the 
more rapidly varying ones to the retarded CEM held of (3.4). 
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\\f DEGREES 

Figure 3: The transverse electric field of a uniformly moving electric charge, scaled by 
the factor r'^/Ze, as a function of angular position i/j, for different rclativistic velocities 
(3 = v/c. The slowly varying sine curves correspond to the RCED formula (3.3), the more 
rapidly varying ones to the retarded CEM held of (3.4). 
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In the frame Sat the transverse impulse^ is given, according to (3.2) as: 

eET{SN)dt' = Ze^ / -^dt' (3.5) 

-oo J—oo T 

Consideration of the component of the electron velocity v perpendicular to r in the scat- 
tering plane (i.e. in the direction of the unit vector ip, see Fig. la) gives the relation: 

= = -V sin (3.6) 

Multiplying both sides of (3.6) by r and noting that, from the geometry of Fig. la, r sin ~ 
b gives, on rearrangement, 

dt' = (3.7) 
bv 

Using (3.7), (3.5) is converted into an integral over ip: 

Ap{Sm) --^l sin Hi^ = (3.8) 

(3.1) and (3.8) give, for the relation between the scattering angle and the impact param- 
eter: 

eiSM) = ^ (3.9) 
pbv 

in agreement, in the small angle limit, with the exact result (i.e. one valid also for 
vc^c) [21]: 

e Ze^ , , 

tan- = 3.10 

2 mv% ^ ' 

Calculating now the transverse impulse in the frame Sg using the RCED formula (3.3) for 

the electric field gives, instead of (3.5): 

/oo roo sin lb 

eEri, Se, RCED)dt = Ze^ / 1^-dt (3.11) 
-oo J—oo T 

Relativistic time dilatation gives the relation^: ^dt — dt' . Hence the right sides of (3.5) 
and (3.11) are identical, leading to the result: 

2Ze^ 

9{Se, RCED) = 9{Sn) = (3.12) 

Using now (3.4) and (3.7) as well as the relation ^ydt — dt', the transverse impulse imparted 
by the LW electric field in the frame Sg is : 

Ap(S„, LW) = r eMS,. LW)dt = ^ f = ^ (3.13) 

J-oo Ybv Jo {1 - (3^ sm^ ip)^ ibv 

(3.13) gives, for the electron scattering angle: 

e{S,,LW) = ^ (3.14) 
7pov 



^Ap is treated throughout as a positive quantity. 

^Note that t is the proper time r in the initial electron rest frame, so that dt' = jdr is the usual time 
dilatation formula relating time increments in 'stationary' and 'moving' inertial frames. 
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Thus for a given value of b the value of the scattering angle is reduced by a factor I/7 as 
compared to the calculation in Sat. 



In Jackson's book [21] a result consistent with (3.12) was found for the calculation in 
the frame Sg using E{Se, LWy. This is because the calculation effectively used, in (3.13), 
the time t' defined in the frame Sjv instead of the correct time t in the frame Sg. Making the 
replacement dt — > dt' in (3.13) and using (3.7) gives a result, the one quoted by Jackson, 
consistent with (3.12). A footnote in [21] indicates that the author has realised that Ap 
was calculated in the frame Se, whereas the relation for the scattering angle, (3.1), is valid, 
instead, in Scm, where the longitudinal electron momentum is '-fmv. There is then the 
statement: 'The reader may verify that (13.3)' (i.e. (3.9) above)'and (13.4)'(i.e. (3.10) 
above) 'are also correct in the frame in which the electron is at rest by transforming the 
angles from the CM system to the laboratory'. Actually, this is impossible to do, because 
the electron scattering angle is undefined in the electron rest frame. The calculations of 
Ap in the frames Sjv and Sg above using the retarded field (3.4) show, instead, results 
that differ by a factor 7. The calculation of the scattering angle performed in Se then 
gives an angle a factor 7 smaller than the known Rutherford scattering result (2.10). Use 
of Ap calculated with (3.4) in Se then leads to a (1/7)^ dependence of the energy-loss 
formula at relativistic velocities, in contradiction both with the Bethe-Bloch formula and 
experiment [24]. 

That the transverse field in (3.3) is larger, by a factor 7, than the static value in (3.2), 
for any value of ip, is a necessary consequence of the invariance of relativistic transverse 
momentum will now be demonstrated. Since the Lorentz transformation leaves invariant 
transverse spatial intervals: Ay — Ay', relativistic transverse momentum must also be 
conserved: 

dy dy' . , 

Pt = m— = rn—— = m'yv — 7777 v (3.15) 

dr dr 

Considering now the transverse momentum, pt, of the electron in Rutherford scattering 
and associating the frame S with Sat and S' with Se^, then, since v' <^ i;, 7' ~ 1. That is, 
the motion of the electron in Se is non-relativistic. If px is generated by transverse forces 
Ft, Flp in S^v, Se respectively, conservation of transverse momentum gives: 

ApT = FtAI = F!j.At' (3.16) 

so that 

^£ = 1^ = ^-7 (3.17) 

Ft At' i ' ^ ' 

since At — At/7 = At' j^' . Thus conservation of transverse momentum requires the 
transverse electric field in Sg to be a factor 7 larger than that in Sat for all values of i). 
This condition is respected by E{RCED) in (3.3), but not by E{LW) in (3.4), for which 
the condition is satisfied only for = n/2. The forces generated by the retarded field (3.4) 
therefore do not, in general, respect the invariance of relativistic transverse momentum, 
in different inertial frames, that is required by SR. 



^The integration variable was a time rather than the angle ip. 

*This, the opposite of the previous assignment, is done for clarity of presentation. Here it is the rest 
and moving frames of the scattered electron rather than those of the nucleus N that are under discussion. 
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4 Breakdown of Gauss' theorem for the electric field 
flux of a moving charge in relativistic classical elec- 
trodynamics 



The total flux of the electric fleld generated by a moving charge: 

(l)E= [ E-dS (4.1) 

is readily calculated from the formulae (3.3) and (3.4) for different choices of a surface, 
S, that completely surrounds the charge. The results obtained for (i) a spherical surface 
centered on the charge and (ii) a cylindrical surface with axis parallel to the direction of 
motion of the charge, of length 2£ and radius a, with the charge at the axis mid-point, 
are, for the RCED and the retarded LW fields: 



Sphere: 



MRCED) = ^[2^+1] (4.2) 




3 

e{LW) = AnZe (4.3) 



Cylinder: 

(t)E{RCED) = 47rZe ((7- i)cos^o + - 1 (4.4) 

V 7 7/ 

4)e{LW) = AttZc (4.5) 

where tan ^0 = 

The flux of the retarded LW fleld is independent of the velocity of the charge and of 
the form of the surface S, but, as can be seen in Fig.3, the electric fleld strength peaks 
strongly near ■0 = 7r/2 for large velocities /3 ~ 1. The number of fleld lines is conserved 
but their spatial distribution changes with the velocity of the charge. This change in the 
distribution of the fleld lines is often interpreted as due to relativistic 'length contraction' 
[25]. However, the flux of the RCED force flelds do not, as is clear from (4.2) and (4.4), 
respect Gauss' theorem, or the electrostatic Maxwell equation, when the source charge is in 
motion. The flux depends both on the velocity of charge and the geometry of the surface 
S. For large values of 7 the flux is proportional to 7, but with different multiplicative 
constants for spherical or cylindrical surfaces. In RCED, both the tranverse fleld and the 
total flux are proportional to 7 at large 7, whereas in the same limit the transverse LW 
fleld is proportional to 7 only for ^/^ = 7r/2, while the total flux is velocity independent. 
There should be nothing shocking in these properties of the flelds. As discussed in [4], the 
force flelds are only mathematical abstractions introduced in order to write inter-particle 
forces in a compact manner. The physical reality resides in the particles and the forces 
between them, not in the flelds. All equations of motion may be written down and solved 
without the necesssity to introduce either the concept of a force or that of a fleld. Just 
this approach is taken in the discussion of Keplerian orbits in Section 6 below. The total 
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flux of the fleld, though readily calculable, has no independent physical meaning and it 
is irrelevant, for the dynamical description of the motion of the charges, whether the flux 
respects, of not, Gauss' theorem. However the forces derived from the fields must respect 
on the one hand, the constraints of SR and on other be consistent, in the appropraiate 
limit, with the known and tested predictions of QED. As the several examples considered 
in the present paper show, this is the case for the instantaneous RCED fields, but not for 
the retarded LW ones. 



5 Resolution of Jackson's torque paradox 



Another example of the incompatibility of the LW retarded fields with relativistic in- 
variance is provided by the 'torque paradox' recently pointed out by Jackson [26]. Replace 
in Fig. lb the moving nucleus by an electron at rest. The electric field at e is then a static 
Coulomb field parallel to r. This is the only force-generating field in the frame Sg and 
gives no torque. Now, on boosting from Sg into Sat, both electrons move with velocity v 
parallel to the x-axis. According to the LW formula (3.4), the electric field is still parallel 
to r and therefore generates no torque. However, in Sjv, the electric current constituted 
by the lower electron produces at e a magnetic field of magnitude: 

B{LW) = /3 X E{LW) (5.1) 

The corresponding Lorentz force on e is parallel to Oy and gives a non- vanishing torque, 
f: 

f=-efx [/3x{f3x E{LW))] = -ke/3^ simp cos^Ij\E{LW)\ (5.2) 

where is a unit vector parallel to the 2;-axis^. Since this torque does not exist in the 
inertial frame Se, related to Sjv by Lorentz transformation, there is a manifest breakdown 
of the relativity principle. This is Jackson's torque paradox. A similar problem was 
previously pointed out [27] in connection with the interpretation of Trouton and Noble 
experiment [28]. 

The RCED fields at e in Sat, on replacing N by an electron moving with velocity v 
along Ox' are, using (2.9) and (2.10), and dropping the primes on the unit vectors in Sn- 

E{RCED) = -^(i^^ + 7jsinV') (5.3) 

B{RCED) = -—7/3 sin V' (5.4) 

The corresponding Lorentz force, F, on e is: 

F{RCED) = -e[E{RCED)+/3xB{RCED)] 
icosip 



+ 7jsin'0 + P(i X A;) (7/3 sin ■0) 



7 

+ J7(l - (3^)sm'ijj 
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^This formula is identical to Eqn(7) of [26] on making the replacements 9 ^ ij), (3q ^ (3 and t = in 
the latter. 
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— - [z cos ip + jsmip] = — -r 



(5.5) 



Thus the Lorentz force in S^v is radial and, consistent with relativistic invariance, there 
is no torque in this frame. 



6 Circular Kepler ian orbits in Relativistic Classical 
Electrodynamics 



In [4] equations of motion, of two mutually interacting charged objects Oi and O2, 
were derived. They may be most simply written in the 'fieldless' and 'forceless' form: 



dt 
dt 



mi 



1712 



d{-fiVi) 
dt 

d{'^2V2) 

dt 



c 



j^f+Pi X (j^ X r) 



]_dj2_ r^ {r-P2) 
cr dt 



02 

c 



j^f+ (52 X (ji X r) ^ 1 dji _ 7^ (f ■ A) 

cr dt 



(6.1) 



(6.2) 



where mi and m2 are the masses of the objects in the absence of mutual interaction, and 
qi and ^2 their electric charges. The vector r = rl2 = i\ — T2 gives the spatial separation 
of the objects in their centre-of-mass frame. The 4-vector current ji of the object Oj is 
defined as ji = QiUi, where Ui is the 4- vector velocity Ui = [c-ffjCjifSi). In the following, 
the motion of the objects in their centre-of-mass frame is considered, and the origin of 
the vectors fi and r2 is the common centre of energy, O, of the objects. 



Because of the relation connecting the velocity, current and energy-momentum 4- 
vectors: 

i n 

(6.3) 



M = - = — 
q m 



the differential equations (6.1) and (6.2) are coupled via the dj/dt terms on their right 
sides. Solving these equations for the time derivatives of the relativistic velocities of the 
two objects gives the equations: 



d{-fiVi) 
dt 



1 



mi 



[1- 



X 



{72[r + /3ix(/32Xr)-/32(f-/32)] 
+^^[f + P2 X (/3i X f) - /3i(f • /3i)] 



(6.4) 



m2C^r ' 



d{'^2V2) 

dt 



-qiq2 



X 



{7i[r + /32x(/3ixr)-/3i(f-/3i)] 
+^^\r + A X 02 X r) - 02{f- P2)] 



(6.5) 



The differential geometry of uniform circular motion is illustrated in Fig.4. The vectors 
Vi or V2 are perpendicular to r at all times. Also \vi\, \v2\, ri and r2 are constant at all 
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times. It follows from the geometry of Fig. 4. that, to first order in 5r,5v : 

50 _ 1 \6vi\ _ 1 \6v2\ _ 1 l^fil _vi_ 1 \Sr2\ _V2 , . 

5t vi 5t V2 5t ri 5t ri r2 5t r2 

For the case of uniform circular motion, qi and q2 must have opposite signs and (6.4) and 
(6.5) simplify to: 



71 



72- 



dvi 
It 

dV2 

It 



-kil|g2|(i + /5i/52) 



mi 


1 


(<Jl'?2)^ 








{1 + I3il32) 


1712 


1 - 


(9l«2)2 





72 



mimni 



7i 



m2C^r 

|gi||g2|72 

m2C^r 



r 



These equations may be combined to yield the relation: 

r.*dvi dV2 _ 

where 



7imi 



^ _ I'?ilk2l7i 

72^2 



7i 



ki||g2|72 



(6.7) 
(6.8) 

(6.9) 

(6.10) 
(6.11) 



Since, for uniform circular motion, £1 and 82 are constant, (6.9) may also be written as: 

j^[£lvi + £;v2] = J^\P\+ p1] = (6.12) 

or, equivalently, 

(6.13) 



dt 



dt 



where p* = S^Vi, i = 1, 2. This is just the expression of Newton's third law for the two 
mutually interacting objects. 

Equations (6.10) and (6.11) then suggest that the electromagnetic interaction between 
the objects modifies their masses as a function of their velocites and spatial separation 
according to the relations: 

mi 



rrir, 



^ _ Igillg2l7i 

1712 



71 



|gl||92|72 



(6.14) 
(6.15) 



Thus, when proper account is taken of the modification of the masses of the objects, due 
to their mutual electromagnetic interaction, there is not, in the case considered here, any 
breakdown of Newton's third law. 



Combining (6.6) and (6.9) 



\dvi\/dt _ 8* 
\dv2\/dt ~ 

14 



Vl 
V2 



H 

r2 



(6.16) 




Figure 4: Differential geometry of tlie position and velocity vectors of two objects, Oi, 
{i = 1,2) executing uniform circular motion about their common center-of-energy O. As 
Sfi and Svi ^ , Sfi±fi, Svi±Vi, \ fi\ — const, and \vi\ — const. 
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or 



Sir, 



(6.17) 
(6.18) 



Equation(6.17), which may be written as \pW — is consistent with momentum conser- 
vation in the overall centre-of-mass frame: p\ = —pi, while (6.18) states that the origins 
of fi and r2 are at the centre of energy of the two interacting objects. 

The period of rotation, r, of either object in its circular Keplerian orbit, is readily 
derived from (6.7) or (6.8) and the geometrical relation (6.6). Considering the object Oi: 



71- 



\dvi 

IF 







k2|(l+/3l/32) 




(5192)2 





72 - 



at 



27rri 27r 



(6.19) 



Using the relations: 



r — n + r2 



ri(l + -) 
ri 



ri 



^^2 



riSl 



where S* is the 'reduced energy': 



(6.20) 



(6.21) 



SI + SI 

to eliminate ri in favour of r on the right side of (6.19) yields a formula relating the 
rotation period of the objects to their spatial separation: 



(27r 


fmi'jiS* 


1 - 




^3 






£t\qi\ 


92|(l+/9l/92) 


1911192171 
'2 m2C-^r 



{27ifS* 


"1 (9192)2 


y-«3 


-L AO 




?ll 


\q2\{l + W 



(6.22) 



where, in the last member of (6.22), (6.10) has been used. This is the relativistic gen- 
eralisation of Kepler's third law. The usual classical result is recovered in the limit 
c — * 00. Note that no kincmatial approximation whatever has been made in the deriva- 
tion. Equation (6.22) gives the exact result for the period of a circular orbit to all orders 
in \qi\\q2\/{mic^r), \qi\\q2\/{m2C^r), and /^a- 



7 The impossibility of uniform circular motion under 
retarded Lieard-Wiechert fields 

Several years before the advent of special relativity the effect on electric and magnetic 
fields, propagated at the speed of light, was considered by Lienard and Wiechert [22]. A 
good explanation of the physics arguments presented by these authors may be found in 
Reference [8]. In this case the values of the fields at a given space-time point are given by 
the sum of the contributions from all sources that lie on the backward three-dimensional 
light-cone of the space-time point considered. Neglecting relativistic effects (as of course 
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Lienard and Wiechert were obliged to do) the instantaneous classical scalar and vector 
potentials 0(^) = q/r and A{X) — qv/r are found to be modified according to the 
relations [8]: 




(7.1) 



(7.2) 



where the quantities in large curly brackets are evaluated at the retarded time: tret = 
t — r/c at which the source charge is sampled by a signal moving along the backward 
light-cone of the space-time point X = {ct;x). The electric and magnetic fields derived 
from these potentials are given by the expressions [29]: 



acceleration respectively of the source charge. 

The possibility of uniform circular motion of two objects with equal and opposite elec- 
tric charges and different masses, interacting mutually via the fields (7.3)-(7.5) above is 
now investigated. As in the previous section, the geometrical constraints of uniform circu- 
lar motion are imposed and a solution of the equations of motion of the objects consistent 
with these constraints is sought. Various geometrical parameters used to describe the sys- 
tem are defined in Fig. 5. The objects Oi and O2 of masses, in the absence of interaction, 
nil and m2 and electric charges q and —q are assumed to perform uniform circular motion 
about their centre of energy O. The radii of the circular orbits of Oi and O2 are ri and r2. 
Cartesian axes Ox, Oy and Oz specified by unit vectors i, j and k are chosen with y-axis 
parallel to r2 at some instant and x-axis in the plane of the orbits. The resultant force on 
Oi due to the retarded fields of O2 is now calculated and examined for consistency with 
uniform circular motion. According to Eqn(7.3) ^^^^(Xi) = -E^gj (1) points to the 'present 
position' (assuming uniform motion) of the source charge. Calculating the retarded fields 
of O2 at Oi when the latter is at position A in Fig. 5 at t = 0, implies that the fields are 
given by the formulae (7.3)-(7.5) when O2 was at position B at time t = —r/c. E^^I{1) 
points towards the point C on the tangent to the orbit of O2 at B. The distance BC is 
r(32^^- The calculation of the retarded fields from the Eqns(7.3)-(7.5) and the geometry 

^°The object O2 does not actually follow this trajectory under its assumed circular orbit. However, 
causality implies that only the velocity of the source at the instant it is sampled on the backward light- 
cone of the point at which the field is evaluated is relevant. Any change in the motion (velocity or 
acceleration) of the source at later times therefore cannot affect the retarded field. The retarded field 
i?^|((l) at A is therefore the same whether O2 follows a circular trajectory BE or the straight one BC. 




(7.3) 



(7.4) 



(7.5) 



E^^l{X) and E^^^{X) give the contributions to the electric field due to the velocity and 
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A 



Figure 5: Geometry of uniform circular motion of two oppositely charged objects, Oi and 
O2, about their common centre of energy O. The forces on Oi at A generated by the 
retarded fields (7. 3)- (7. 5) of O2 are calculated. The force fields are produced by O2 at 
the position B. Various geometrical parameters are defined. The directions of the electric 
fields at A associated with the velocity and acceleration of O2 are shown in the inset 
figure. 
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of Fig. 5 is presented in the Appendix. The following results are obtained: 



4"e1(l) = 



iq/3l , ig[l-/52(2/5i + /52/2)] 



3r2 



+ 



r2(ri + r2 



in + r2)^ 
+ 0(/32') 



+ 0(/32') 



+ 0{(3t) 



(ri + r2) 

According to the Lorentz force equation, the total force acting on Oi is: 



(7.6) 
(7.7) 
(7.8) 



Fret{l) 



r2 



A + P2^(32_ 

ri + r2 3r2 



+ 



ri + r2 



l-/32(2/3i + /32/2) 
n + r2 



(7.9) 



By symmetry, the force on O2 due to the retarded fields of Oi is: 



Freti'2) 



I3i+I32 
ri + r2 



3ri 



ri + r2 



/3i(2/32 + A/2) 
ri + r2 



+ 0{Pt) (7.10) 



Inspection of (7.9) and (7.10) demonstrates the impossiblity of uniform circular motion 
under the action of retarded Lienard - Wiechert fields. There is a manifest breakdown 
of Newton's third law, that cannot, due to the non-factorisable force components on the 
right sides of (7.9) and (7.10), be restored, as in the discussion of the previous section, by 
the introduction of different effective masses for the interacting objects. Also, since the 
resultant force is non-central, an unbalanced torque acts on the system. The unbalanced 
force is: 



+ 



' Pl + f32 

ri + r2 

PI' 



+ 



2(ri +r2)2 



(7.11) 



The clockwise (accelerating) torque is of magnitude: 



2(/3i + /32) I 1 //32 ^ 
ri + r2 3 I r2 n . 



(7.12) 



In [19] Eddington discussed the case of the (almost circular) orbits of Jupiter and the Sun 
under the action of retarded gravitational fields. A simple geometrical argument suggests 
that Jupiter and the Sun would then be subject to accelerating torques. Eddington 
then considered the analogous electromagnetic problem, and citing the Heaviside formula 
(3.4) for the electric field of a moving charge, stated that the direction of the retarded 
force is 'very nearly' in the direction of the centre of mass in this case. It was then 
concluded that there is no unbalanced torque problem for the Sun- Jupiter system, and 
that it may be generally assumed that the gravitational force propagates at the speed 
of light. The calculation presented above shows this conclusion to be incorrect. As the 
formulae (6.1) and (6.2) follow from only the inverse- square force law, SR and Hamilton's 
Principle, similar formula, with the relacement qiq2 — > —Gmim2, are expected to hold 
in Newtonian gravitation, when special relativistic corrections are included. In this case 
also, stable circular orbits cannot be obtained with retarded gravitational fields. 
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A calculation considering the stability of circular orbits for the case of gravitational 
forces has been previously published [30]. As for the electromagnetic forces just consid- 
ered, and for similar geometrical reasons, an accelerating torque was found in the case of 
retarded gravitational fields. By comparing the predicted rate of change of their periods 
with precision data on the rate of change of the observed periods of binary pulsars, a 
lower hmit of 2800c for the speed of propagation of gravitational force fields was set. The 
conclusions of Ref. [30] were questioned in Refs. [31, 32]. However, the argumentation of 
these authors was based on an analogy between retarded gravitational and electric force 
fields. They concluded that since the retarded LW electric field of Eqn(3.4) is radial at 
the 'present time' there is no torque. As discussed above, a similar remark was made 
by Eddington [19]. This conclusion is refuted by the calculation presented above. The 
present-time retarded electric field is radial only for a source charge in straight-line mo- 
tion. As the above calculation and the analogous gravitational one of Ref. [30] show, the 
force is no longer radial when the circular geometry of the orbits is properly taken into 
account. Carlip's argument in Ref. [31] has also been rebutted in Ref. [33]. 



8 Summary 

The present paper is a sequel to [4] where inter-charge forces in classical electrody- 
namics were derived from Coulomb's law, SR and Hamilton's Principle. An essential 
aspect of the derivation is the instantaneous nature of the electromagnetic forces, demon- 
strated in [4] to be a prediction of QED. In standard text books on classical electromag- 
netism [8, 9, 10] no clear distinction is made between fields describing mechanical forces 
(the effect, in QED, of the exchange of space-like virtual photons) and fields giving a 
classical description of radiative processes (in QED, the creation or destruction of real, 
on-shell, photons). The importance of such a distinction, the former fields being instan- 
taneous, the latter propagating at the speed of light, was pointed out many years ago 
by Chubykalo and Smirnov-Rucda [34]. This conjecture is fully supported by the work 
presented in [4] and the present paper. 

The distinction between the standard text-book formulae, denoted in the present paper 
by the labels CEM or LW (for Lienard-Wicchcrt) and those (denoted by RCED) derived 
in [4] or Section 2 above, is most succinctly expressed in terms of the electromagnetic 
potentials generated by uniformly moving electric charges in the different theories. In 
CEM, no distinction is made between fields describing mechanical forces and radiation 
processes. In both cases the electromagnetic potential is assumed to be of the form first 
derived by Lienard and Wiechert [22]: 

i(C.M) - (8.) 

where all quantities in the large curly brackets are evaluated at the retarded time tret = 
t — r/c. In RCED there are two physically distinct types of fields. The first, describing 
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mechanical inter-charge forces are derived, using (2.1) and (2.2), from the potentials: 



qu 







A'f^iRCED) 
AforiRCED) 



r cr 
q^v qu 



cr 



cr 



(8.3) 



(8.4) 



where u = (7c : ^v) is the 4- vector velocity of the moving charge. These potentials as well 
as the associated fields and forces are instantaneous. The retarded potentials of RCED, 
describing radiation, arc the same as (8.3) and (8.4) above, except for a retarded time 
argument. The different physical interpretations of the fields derived from A for and the 
retarded potential are discussed in [4]. See also [34]. It is perhaps surprising, 100 years 
after the advent of SR, that the formulae (8.1) and (8.2) are still the generally accepted 
ones for the electromagnetic potential of a moving charge. Indeed A^ and A in these 
equations are manifestly not components of a 4-vector, and so it is not to be expected 
that the fields and forces derived from them will respect SR. It has been demonstrated in 
the present paper that indeed they do not, whereas the RCED fields and forces do. 

In fact since (8.1) and (8.2) were derived by Lienard and Wiechert some seven years 
before Einstein's first SR paper, it would be little short of miraculous if they constituted 
a 4-vector. A similar remark applies to the consistency with SR of the 'present time' 
formula (3.4) for the retarded electric field of a uniformly moving charge. This may be 
derived from (8.1) and (8.2), but was first obtained by Heaviside [23] in 1888. It may 
be noted, on the other hand, that a velocity-dependent scalar potential, consistent with 
(8.3) up to 0(/9^) was proposed as early as 1861 by Riemann [35]^^. 

The mechanical forces in CEM derived from (8.1) and (8.2) have been compared above 
with those of RCED derived from (8.3) and (8.4). In Section 3 it is shown that use of 
the CEM formula (3.4) to calculate small-angle Rutherford scattering gives a scattering 
angle smaller by a factor I/7 than the well-known standard result, that is recovered by 
use of the RCED formula (3.3). It is also shown in Section 3 that that conservation 
of relativistic transverse momentum requires that transverse forces in different inertial 
frames scale according to I/7 (see Eqn(3.17)). This condition is respected by the RCED 
fields but not by the CEM ones. In Section 4 it is shown that the flux of the electric 
field of a moving charge in RCED does not (unlike in CEM where the electric flux is the 
same for a charge at rest or in motion) respect Gauss' Law. The flux depends both on 
the velocity of the charge and the geometry of the bounding surface. The electric field 
is a mathematical abstraction useful to simplify dynamical formulae expressed in terms 
of forces, but has, unlike charges and their motion, no objective existence. In Section 5, 
Jackson's 'torque paradox' [26] for two interacting charges, where electric and magnetic 
fields are derived from the CEM potentials (8.1) and (8.2) is discussed. A non-vanishing 
torque is predicted in one inertial frame, a vanishing one in another, in contradiction with 
the special relativity principle. It is shown that use of RCED fields gives a vanishing 

^^The same Bernhard Riemann (1826-1866) is better known for his hypothesis concerning the distri- 
bution of prime numbers, for pioneering applications of topology to the theory of complex functions and, 
with Lobatchevsky, the invention of non-Euclidean geometry. Together with Weber and Clausius, he 
attempted to construct an action-at-a-distance theory of electromagnetic forces in which the fundamen- 
tal physical objects were charged particles, not fields. Sec the second reference in [35]. This is also the 
approach adopted in [4] by the present author. Such theories were generally superseded in the second 
half of the 19th century by Maxwell's electromagnetic field theory. 
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torque in all inertial frames, consistent with SR. In Section 6 it is demonstrated that 
stable, circular, Keplerian orbits are possible for two oppositely-charged objects using 
exact RCED dynamical formulae. An exact relativistic generalisation of Kepler's third 
law relating the orbit period to the separation of the charges is derived for this case. In 
Section 7 the impossiblity of stable circular motion under forces derived from retarded 
CEM fields is demonstrated. Such forces do not respect Newton's third law and, as 
already conjectured during the 19th Century for the analogous gravitational case [19] and 
recently demonstrated by explicit calculation [30] , generate an unbalanced torque. Stable 
orbital motions of electrons in atoms or of planets around the Sun are only possible if (as 
Newton assumed for the case of gravity) the inverse-square-law electrical or gravitational 
forces act instantaneously. 
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Appendix 

Denoting by 5t the time of passage of O2 from B to E (Fig. 5) then: 
St^'-^^ (Al) 

C V2 

Prom the geometry of Fig. 5: 

r — ri cos 01 + r2 cos 02 (^2) 
ri sin 01 — r2 sin 02 (^3) 

= 01 + 02 (^4) 

riV2 = r2'i;i (^5) 
It follows from (A1-A5) that: = /3i + /32, 0i = /32+O(0?) and 02 = /3i+O(0i) 
Making the small angle approximations: 

03 , 02 

sin ~ — — cos ~= 1 — — 
o ^ 

equations (A2)-(A5) may be combined to give: 

r ~ (ri + r2)(l-^) 
2r2 

= (ri + r2)(l-^) + 0(/3|) (^6) 



r-(32-r = r(l + /32sin02) 
~ r(l + /3i/32) 



(n + r2)(l + ^) + 0(/3,3) (A7) 



also 



= f — /^r = r'(zsinQ; — jcosa) (A8) 
Substituting (A6)-(A8) into (7.3) and noting that 52 = —1 and r' = rsec(0i — a) gives: 



?veiM\ g(l -/3|)(-zsina + jcosa) 



r2(l + /?i/32)3cos(0i-Q;) 

q{l — /3|)(1 + ^)(— zsino; + jcosct) 

(ri + r2)2(l + 2/?i/52) 



g[-m + i{l-/j2(2A + /32/2)}] 4 .4^ 

The last member follows since, as shown below, a is of order 
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Prom the geometry of Fig. 5: 

CD cos ( 



tan a 



Ti + r2 sec (j) — CD sin 

(r2 tan 6 — ri^2) cos (6 
r\ + r2 sec — (r2 tan — rli^) sin 

r2(sin0 — 0COS0) 
ri + r2 sec — r2(tan — 0) sin 

where (Al) has been used. Noting that 

6^ 6^ 

sin0-0cos0 = 0-^-0+^ + ... = ^ + 0(05) 

D Z D 

tan0sin0 - 0sin0 = ^'^ ~ - 02 + ^ + ... = ^ + 0(0^) 

^ 2" 



sec0 = 1 + — + ... 



(AlO) gives the result: 



q; ~ arctan 



r2 



rx^r-i) 3 



A3 



r2 / 3 



^ + (All) 



Substitution of a. given by (All) into (A9) gives Eqn(7.6) of the text. 

Since the modulus of the velocity vector of an object in uniform circular motion is 

— * 

constant, the increment b^i for O2 is anti-parallel to the radius vector r2. It then follows 
from Fig.5 and the geometrical conditions (6.6) that: 

fx [(f - /fsr) X dMdt\ = !I^^!l^[£cos0i + isin0i] + Oi^^) iAVl) 

Subsituting (A12) and (A7) into (7.4) gives: 

q(5l {Pi +P2)[i cos 01 + j sin 0i] 



^."eT(l) 



rr2(l + cos(0i - a) 

2 



?/3|(l + |)(/3i + /32)[£(l-|)+i/?2] 



(ri + r2)r2(l + /3i/32)=^(l-^) 
which yields Eqn(7.7) of the text on retaining only terms up to 



{A13) 
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